Understanding Quantum Entanglement: Qubits, Rebits and the 

Quaternionic Approach 



J. Batle^*, A. R. Plastino^-^-s, M. Casas^ and A. Plastino^'^ 
^ Departament de Fisica, Universitat de les Illes Balears, 07071 Palma de Mallorca, Spain 
'^National University La Plata, C.C. 727, 1900 La Plata, Argentina 
\ ^Argentina's National Research Council (CONICET) 

;h : (February 1, 2008) 

Abstract 

> 
O 

O . 

^ . It has been recently pointed out by Caves, Fuchs, and Rungta [1] that 

. real quantum mechanics (that is, quantum mechanics defined over real vector 



X 



spaces [2-5]) provides an interesting foil theory whose study may shed some 



. light on just which particular aspects of quantum entanglement are unique to 

^ : 

^ . standard quantum theory, and which ones are more generic over other physical 

theories endowed with this phenomenon. Following this work, we discuss 
some entanglement properties of two-rebits systems, making a comparison 
with the basic properties of two-qubits systems, i.e., the ones described by 
standard complex quantum mechanics. We also discuss the use of quaternionic 
quantum mechanics as applied to the phenomenon of entanglement. 
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I. INTRODUCTION 



The phenomenon of entanglement is one of the most characteristic non-classical features 
exhibited by quantum systems [6]. Quantum entanglement is the basic resource of sev- 
eral quantum processes as superdense coding [7], quantum teleportation [8] and quantum 
computation [6,9] studied by quantum information theory [10-13]. A state of a composite 
quantum system constituted by subsystems A and B is called "entangled" if it cannot be 
represented as a mixture of factorizable pure states. 

The simplest systems exhibiting the phenomenon of entanglement are two-qubits sys- 
tems. For these systems an explicit expression of the entanglement of formation, which 
is a natural quantitative measure of entanglement [14], has been found by Wootters [15]. 
The correspondent space of mixed two-qubits states in standard quantum mechanics is 15- 
dimensional. The amount of entanglement in this space has been established by Wootters 
[16-19]. 

For quantum mechanics defined over real vector spaces the simplest composite systems 
are two-rebits systems. An explicit expression for the entanglement of formation of arbitrary 
states of two-rebits has been obtained by Caves, Fuchs and Rungta [1]. Pure states of rebits- 
systems are described by normalized vectors in a two dimensional real vector space. The 
correspondent space of mixed two-rebits states is 9-dimensional (vis-a-vis 15 for 2 qubits). 

The aim of the present work is to explore numerically, as well as conceptually, the 
entanglement properties of two-rebits systems [20], as compared to the usual two-qubits 
ones, so as to detect the differences between the two types of system. We pay particular 
attention to the distribution of states according to their degree of entanglement. We obtain, 
analytically the probability densities for finding pure states states with a given amount of 
entanglement E (or with fixed values of the so called concurrence squared C^). The same 
is done for mixed states in numerical fashion. 

The paper is organized as follows. In section II we review the general properties of two- 
qubits systems. Several quantities related to the entanglement of formation are investigated 
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for both rebits and qubits in section III, with an emphasis on the differences and simihtudes 
between both formahsms. In section IV we discuss how to follow a quaternionic approach 
to quantum entanglement. Finally, some conclusions are drawn in section V. 

II. ENTANGLEMENT OF TWO-QUBITS SYSTEMS 

As already stressed in the introduction, the two-qubits system is the simplest quantum 
mechanical system that exhibits the sort of "extra correlation" one calls entanglement be- 
tween two parties. A representation of a qubit is given by the Bloch sphere. The poles 
correspond to classical bits |0), |1) while any point on this unit sphere, given by a pair of an- 
gles (p, ip, represents a general qubit cos0|O) -|-e*'^sin0|l). Qubits constitute the essential new 
elements in the switch from Classical to Quantum Information Theory and entanglement 
the basic ingredient of many striking processes now under active investigation. 

Here we use the concurrence C[p] and the entanglement of formation E\p] as quantita- 
tive measures of the amount of entanglement. Wootters' formula for the entanglement of 
formation of a two-qubits state p reads [15] 

where 

h{x) = — a;log2a; — (1 — a;) log2(l — a;), (2) 

and the concurrence C is given by C = 'max{0, Ai — A2 — A3 — A4), with Aj, (i = 1, . . . , 4) 
being the square roots, in decreasing order, of the eigenvalues of the matrix pp, with 

p = {ay(»ay)p*{ay(^ay). (3) 

The above expression has to be evaluated by recourse to the matrix elements of p computed 
with respect to the conventional product basis. 

In this work we have generated all pure and mixed states of a two-qubits system according 
to the measure defined in [16,17,21]. 
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The distribution of two-qubit states according to their degree of mixture can be obtained 
analytically adopting a geometric representation [22] for the eigenvalues of p as a regular 
tetrahedron Ta of side length 1, in 7?.^, centered at the origin. There exists a mapping 
connecting these eigenvalues (pi, . . . ,P4) with the points r of the tetrahedron, so that we 
can relate the participation ratio R{p) = 1/Tr(p^) to the radius r =| r | of a sphere embedded 
within the tetrahedron Ta according to 

r'--l + lf.Pl (4) 

Thus, the states with a given degree of mixture he on the surface of a sphere of radius r 
concentric with the tetrahedron Ta- For a completely description, see [22]. It is then easy to 
find out just how our states are so distributed, according to the degree of mixture R. This 
distribution is depicted in Figure 1. The dualistic nature of the concepts of entanglement 
and mixedness becomes apparent. For two-qubits systems one would expect that the states 
will tend to become less entangled as the degree of mixture increases. In fact, for i? > 3 
(vertical line in Figure 1) all states are separable. It is important to stress here that this 
curve is the same no matter what quantities (rebits, qubits, or even quaterbits) are involved, 
although their respective entanglement properties are different. 

III. ENTANGLEMENT PROPERTIES OF TWO-REBITS SYSTEMS. A 

COMPARISON 

In the space of real quantum mechanics wc can represent rebits on the Bloch sphere. The 
poles correspond to classical bits 1 0) , 1 1) , but the sphere reduces itself now to a maximum unit 
circle, described by just one parameter 0. We have cos0|O) +e''^sin0|l) cos0|O)+sin0|l). 
Entanglement can also be described in such a context with suitable modifications. Caves, 
Fuchs, and Rungta's (CFR) formula for the entanglement of formation of a two-rebits state 
p is given by (1), (2) using the concurrence C[p] — | tr(T) | = | tr{pay (8) ay) |, which has to 
be evaluated using the matrix elements of p computed with respect to the product basis, 
=1 iJ = 0,1. 



For a two-rebits state the entanglement of formation is completely determined by the 
expectation value of one single observable, namely, ay ® (7y. On the contrary, it has been 
recently proved that there is no observable (not even for pure states) whose sole expectation 
value will yield enough information to as to determine the entanglement of a two-qubits state 
[23] . As shown in [20] , there are mixed states of two rebits with maximum entanglement (that 
is, with = 1) within the range 1 < R < 2. This is clearly in contrast to what happens 
with two-qubits states, because only pure states {R — 1) have maximum entanglement. 

The measure of mixedness R introduced above can be used in the description of the 
entanglement phenomenon for two-rebits systems. Also, in order to explore numerically the 
properties of arbitrary two-rebits states, it is necessary to introduce an appropriate measure 
11 on the space Sr of general two-rebits states. Such a measure is needed to compute 
volumes within the space Sr, as well as to determine what is to be understood by a uniform 
distribution of states on Sr [16,17]. 

An arbitrary (pure and mixed) state p of a (real) quantum system described by an 
A'"-dimensional real Hilbert space can always be expressed as the product of three matrices, 

p = RDliXiWR". (5) 

Here R is a.n N x N orthogonal matrix and D[{Ai}] is an iV x diagonal matrix whose 
diagonal elements are {Ai, . . . , A^v}, with < Aj < 1, and — The group of orthogonal 
matrices 0{N) is endowed with a unique, uniform measure v [21]. On the other hand, the 
simplex A, consisting of all the real A^-uples {Ai, . . . , Ajv} appearing in (5), is a subset of a 
{N — 1) -dimensional hyperplane of TZ^ . Consequently, the standard normalized Lebesgue 
measure jCn-i on TZ^~^ provides a natural measure for A. The aforementioned measures 
on 0{N) and A lead then to a natural measure p = u x C^^i on the set Sr of all the states 
of our (real) quantum system. 

Clearly, our system will have N — A again, and in our numerical computations we are 
going to randomly generate states of a two-rebits system according to the measure p. 

The relationship between the amount of entanglement and the purity of quantum states 
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of composite systems has been recently discussed in the hterature [16-19]. The amount 
of entanglement and the purity of quantum states of composite systems exhibit a dualistic 
relationship. As the degree of mixture increases, quantum states tend to have a smaller 
amount of entanglement. In the case of two-qubits systems, states with a large enough degree 
of mixture are always separable [16]. To study the relationship between entanglement and 
mixture in real quantum mechanics, we compute numerically the probability P{E) of finding 
a two-rebits state endowed with an amount of entanglement E. In Figure 2 we compare 
(i) the distribution associated with two-rebits states with (ii) the one, associated with two- 
qubits states, recently obtained by Zyczkowski et al. [16,24]. Fig. 2a depicts the probability 
P{E) of finding two-qubits states endowed with a given entanglement E (as computed with 
Wootters' expression). In a similar way. Fig. 2b exhibits a plot of the probability P{E) 
of finding two-rebits states endowed with a given entanglement E (as computed with the 
CFR formula). Comparing Figs 2a and 2b we find that the distributions P{E) describing 
arbitrary states (that is, both pure and mixed states) exhibit the same qualitative shape for 
both two-qubits and two-rebits states: in the two cases the distribution P{E) is a decreasing 
function of E. 

The distribution P{E) or P{C^) for pure two-rebits states can be obtained analytically. 
Let us write a pure two-rebits state in the form 

I *) = E Q I <f>i), (6) 

i=l 

where 

= 1, Q e 7^. (7) 

i=l 

The states (| i = 1, . . . , 4) are the eigenstates of the operator Oy^Gy. The four real 
numbers q constitute the coordinates of a point lying on the three dimensional unitary hyper- 
sphere S's (which is embedded in V}). We now introduce on ^3 three angular coordinates, 
01, 02, and 6, defined by 

Ci = cos ^ cos 01, 
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C2 = COS ^ sin 01, 
C3 — sin ^008 02, 

71" 

C4 = sin 61 sin 02, <0 < -, < 0i,02 < 27r. (8) 

In terms of the above angular coordinates, the concurrence of the pure state | ^) is given by 

C =|(cTj^(g)crj^)| = |cos26'| . (9) 

Using (8) and (9) one deduces that the probabihty density P(C^) of finding a pure 
two-rebits state with a squared concurrence is given by 

P(C2) = (10) 



The distribution is to be compared with the one obtained for pure states of two-qubits 
systems, which is (analytically) found to be [24] 



p(c') = ^vT^. (11) 

Both distributions are compared in Figure 3. Figure 3a depicts the one for qubits, while 
Fig. 3b shows the one for rebits. The distribution remains finite, in the case of qubits, for 
all C^. In the case of rebits it presents a sharp peak at the origin, an then saturates to 
1/2 at = 1. The general conclusion that one draws from Figures 2 and 3 is that the 
curves representing the distributions P(-E) and P(C^) associated with (i) pure states and 
(ii) arbitrary states do not differ, in the case of two-rebits states, as much as they do in the 
case of two-qubits states. 

We can determine analytically which is the maximum entanglement of a two-rebits 
state compatible with a given participation ratio R. Since £^ is a monotonic increasing 
function of the concurrence C, we shall find the maximum value of C compatible with a 
given value of R. In order to solve the ensuing variational problem (and bearing in mind 
that C —\{oy (g) (Ty)| ), let us jirst find the state that extremizes Tr(p^) under the constraints 
associated with a given value of I^Oy ® ay), and the normalization of p. This variational 
problem can be cast in the fashion 



5[lV(p2) + P{ay ay) - aTr(p)] = 0, 



(12) 



where a and (3 are appropriate Lagrange multipliers. 

After some algebra, and expressing the expectation value of {ay ® ay) in terms of the 
parameter one finds that the maximum value of compatible with a given value of R 
is given by 



IV. THE QUATERNIONIC APPROACH TO QUANTUM ENTANGLEMENT 

The quaternionic space 7i constitutes a generalization of the complex space C which, in 
turn, generalizes the real space TZ. Each step of this chain is possible by introducing new 
quantities: i'^ — —1 from TZtoC and j'^ — k'^ — —1 from C to H, with suitable commutation 
laws for the three quantities i,j,k. 

A general quaternion (p and its associated commutation algebra are written in the fol- 
lowing fashion 



with (f)i E TZ and ij — —ji — k, jk — —kj — i, ki — —ik — j. 

This "natural extension" of complex numbers that yields quaternions cannot be gener- 
alized any further. Thus, if we give up the property of commutativity, the most general 
algebra that can be used in quantum mechanics is the quaternionic one [25] . 



The definition of entanglement in Quaternionic Quantum Mechanics (QQM) for pure 
states does not differ from the standard one. Given a pure state \ip) of a composite bipartite 
system, the entanglement is obtained via the von Neumann entropy of the marginal density 




(13) 



(f) ^ (/)Q + i(f)i + j(f)2 + k(f)3, 



(14) 



A. Entanglement for pure states of two "quaterbits" systems 
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matrix associated to subsystem A by tracing over the subsystem B: pA —'^B\i^)AB{i^\ of the 
total density matrix p = \'4')ab{'4'\j or, vice versa, pb —TrA\ip)AB{'ip\- Thus, E{p) = S{pa) = 
S{pb)- 

In the case of quaternions we face a higher dimensionahty and, therefore, we need more 
parameters to describe the state p. Additionally, in using the kets and bras notation of 
Dirac's we must keep in mind the quaternions' non-commutativity rules. For the sake of 
simplicity, let us suppose that a pure state is written as 



\iP)^Ci\0)a\0)b + C2\1)a\1)b; Ci,C2en, |Ci|2 + |C2|' = l. 



(15) 



The statistical matrix p — \iP)ab{'4^\ reads 








Ci Ca 



^ CiC~i CiC^ ^ 





C2CT CsC^ 



(16) 



with p^ — p since C1C2 — C2C1 — C2C1. Entanglement then is only a function of N(CiY — 
Cfo + Cf^ + + Cfr^- Identifying x = Cf^ + Cf^ and y = Cl^ + C13, we plot E{p) in Fig. 
4, together with E[p) for the three Quantum Mechanics' versions we are dealing with here. 



B. Extension to general mixed states 

The full analytical study of entanglement in the framework of Quaternionic Quantum 
Mechanics requires careful consideration of the algebra of states and operators for these 
"hyper- numbers" . If one wishes to discuss how to carry out statistical studies and how 
entanglement-related properties are distributed over the space of all (pure and mixed) states, 
one notices that in this case the dimensionality of the problem for a general 4x4 matrix is 
substantially higher (3-1-4*6=27) than for the complex (3-1-2*6=15) or the real (3-1-1*6=9) 
cases. The ensuing statistical properties become clearly non-trivial and some substantial 
effort is required. 
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Let us merely list here the basic ingredients needed for a complete description of the 
statistical properties of quaternionic states p: i) to build and correctly parametrize the 
unitary transformations defined over the quaternionic Hilbert space, ii) to specify the form 
of the Haar measure for the concomitant space of unitary transformations. Notice that 
the measure on the simplex is exactly the same as for complex or real systems. For mixed 
states p, the distributions associated to quantities that depend only on the eigenvalues of a 
statistical operator do have the same form in the real or the complex cases, for they only 
depend on the simplex. The results of the corresponding numerical study will be published 
elsewhere. 

V. CONCLUSIONS 

We have explored numerically the entanglement properties of two-rebits systems. A 
systematic comparison has been established between many statistical properties of two- 
qubits and two-rebits systems. We paid particular attention to the relationship between 
entanglement and purity in both quantum mechanical frameworks. We have also determined 
numerically the probability densities P{E) of finding (i) pure two-rebits states and (ii) 
arbitrary two-rebits states, endowed with a given amount of entanglement E or concurrence 
squared C^. In particular, we determined analj^ically the maximum possible value of the 
concurrence squared of two-rebits states compatible with a given value of mixedness R. 
As for the probability of finding states with a given amount of entanglement, the difference 
between mixed and pure sates is much larger for qubits than for rebits. Also, we have 
sketched the manner in which the quaternionic formalism could be applied to the study of 
quantum entanglement. 
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FIGURE CAPTIONS 



Fig. 1- Probability (density) distribution for finding a state p witfi a given participation 
ratio R. States (two-qubits) with R>2> are always separable. 

Fig. 2- Probability (density) distributions for finding a state p with a given amount of 
entanglement E. Fig. a corresponds to qubits, while Fig. b to rebits. 

Fig. 3- Probability distributions for finding a state p with a given concurrence squared C^. 
Fig. a corresponds to qubits, while Fig. b to rebits. 

Fig. 4- Plot of the entanglement of pure states in the framework of Real Quantum Mechan- 
ics (RQM), Standard Quantum Mechanics (SQM) and Quaternionic Quantum Mechanics 
(QQM). See text for details. 
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